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ABSTRACT 

Magnetic reconnection plays an essential role in the generation and evolution of astrophysical 
magnetic fields. The best tested and most robust reconnection theory is that of Parker and Sweet. 

1/2 

According to this theory, the reconnection rate scales with magnetic diffusivity as \^ . In the 
interstellar medium, the Parker-Sweet reconnection rate is far too slow to be of interest. Thus, a 
mechanism for fast reconnection seems to be required. We have studied the magnetic merging of 
two oppositely directed flux systems in weakly ionized, but highly conducting, compressible gas. 
In such systems, ambipolar diffusion steepens the magnetic profile, leading to a thin current sheet. 
If the ion pressure is small enough, and the recombination of ions is fast enough, the resulting 
rate of magnetic merging is fast, and independent of Xq. Slow recombination or sufficiently 
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large ion pressure leads to slower merging which scales with An as . We derive a criterion 
for distinguishing these two regimes, and discuss applications to the weakly ionized ISM and to 
protoplanetary accretion disks. 

Subject headings: conduction — diffusion — ISM:magnetic fields — MHD 



1. Outline of the Problem 

The Ohmic diffusion time scale in the interstel- 
lar medium (ISM) is so long compared to the dy- 
namical time scale that at first sight the magnetic 
field would appear to be perfectly frozen to the 
plasma. The ratio of these two time scales, - the 
Lundquist number 5 - is typically 10^^ < S < 
10^^, implying Ohmic diffusion times much longer 
than the age of the universe. 

The weakness of Ohmic diffusion causes seri- 
ous problems for Galactic dynamo theory, which 
requires breaking the frozen-flux condition in or- 
der to convert small scale fields to large scale 
fields. More generally, it is difficult to reconcile the 
frozen-flux condition with the apparent smooth- 
ness of the magnetic field in the ISM and the ap- 
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parently turbulent nature of the interstellar veloc- 
ity field. 

Magnetic reconnection is a "hybrid" process 
which combines resistive and dynamical effects, 
and acts at a rate intermediate between them. 
When the reconnection rate depends on S (usu- 
ally a power-law dependence) , it is said to be slow. 
Fast reconnection, by definition, proceeds at a rate 
independent of 5. It is generally accepted that 
in order for reconnection to cause a substantial 
breakdown of flux freezing in the ISM, it must be 
fast. 

The reconnection theories of Parker (1957) and 
Sweet (1958) have proven durable, and represent 
in some sense the "default" models, to which later 
work is usually compared. In the Parker-Sweet 
models, two opposing magnetic field systems are 
pressed together, forming a current sheet at the 
magnetic null plane at which the field undergoes 
steady state, resistive annihilation. According to 
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this model, the reconnection rate scales as S^/'^. 

Petschek (1964) argued for changing the ge- 
ometrical setup from reconnection in a current 
sheet to an X-point geometry. In this model, 
the reconnection rate is proportional to (In 5)"^. 
However, both numerical (Cowley 1975; Biskamp 
1986; Scholer 1989; Biskamp 1994; Ugai 1995) and 
analytical studies (Kulsrud 1998, 2001) cast doubt 
on whether Petschek's model is realizable, unless 
special conditions are met (Forbes 2001). Another 
solution based on flow geometry was proposed by 
Lazarian & Vishniac (1999), who argued that fast 
reconnection occurs in the presence of turbulence 
at the resistive length scale; see also Chiueh & 
Zweibel (1989). 

The solution to the fast reconnection problem 
could lie in physical processes outside the scope 
of single-fluid magnetohydrodynamics. The high 
electric current density expected in the reconnec- 
tion layer may lead to instabilities which could 
generate small scale electromagnetic or electro- 
static turbulence. Scattering of electrons by such 
turbulence provides so-called anomalous resistiv- 
ity (see Treumann (2001) and references therein) 
and increases the reconnection rate. If the current 
layers are sufficiently thin, electrons and ions may 
decouple (Mandt et al. 1994). As Biskamp et al. 
(1995) argue, this can lead to reconnection at a 
rate independent of the resistivity. 

Dorman & Kulsrud (1995) showed that the re- 
connection layer can be made thin if the gas pres- 
sure is reduced by progressive cooling. They found 
that a resistive region will form within a cooling 
time multiplied by InS, and that the magnetic 
field is dissipated within a few subsequent cooling 
times. However, as they argue, the temperature 
must be reduced by a factor of 5^/* for this to hap- 
pen. Bearing in mind the value of S in the ISM, 
the study demonstrates very clearly how large S 
actually is. 

In this paper, we consider ambipolar diffusion 
(Mestel & Spitzer 1956) as a mechanism to re- 
duce the reconnection layer width. Zweibel (1989) 
investigated the effect of partial ionization on 
tearing instabilities in a sheared magnetic field. 
The reconnection rate is increased by a factor of 
(p/PiV^^j but the S''^/'''-dcpendcncc foimd in the 
theory of Furth et al. (1963) is unchanged. Vish- 
niac & Lazarian (1999) considered Parker-Sweet 
reconnection in weakly ionized gas. They argued 



that although the rate remains proportional to 
5^/^, it can be enhanced by several orders of mag- 
nitude due to recombination of ions in the resis- 
tive layer. Brandenburg & Zweibel (1994) showed 
that ambipolar diffusion steepens the magnetic 
profile in the vicinity of a magnetic null layer, 
and suggested that the resulting high current den- 
sity could set the stage for magnetic reconnec- 
tion. Brandenburg & Zweibel (1995) - hereafter 
BZ95 found that although ion pressure tends to 
broaden the resistive layer and reduce the recon- 
nection rate, a fast reconnection regime exists that 
is accessible under ISM conditions. 

This paper is in some respects a foUowup of 
BZ95. We consider the steady state properties of 
a layer of weakly ionized gas around a field re- 
versal, with the magnetic field fixed at the outer 
boundaries. We show by analytical arguments and 
a set of numerical models, that magnetic merging 
can be either fast or slow. The merging rate de- 
pends on a single parameter combining the ratio 
f3 of ion pressure to magnetic pressure, the Ohmic 
diffusivity An, the ambipolar diffusivity Xad, the 
recombination rate of ions, and the global length 
scale L in the problem. Our result suggests that 
fast merging of magnetic neutral sheets is possible 
in the weakly ionized interstellar gas. 

In §2 we define the problem, quantify the rele- 
vant time scales, and develop an analytical model 
for predicting the reconnection rates. A numerical 
solution of the problem requires a solver which is 
able to handle the disparate time scales; we de- 
scribe the solver in §3. In §4, the numerical re- 
sults are analyzed and compared with the analyt- 
ical theory. We discuss the implications for the 
ISM in §5. Section 6 is a summary and discussion. 

2. Analytical Description 

Before embarking on a mathematical represen- 
tation of the problem, we describe the physical 
setup and identify the critical parameters which 
we would like to determine. 

We begin by summarizing reconnection theory 
of Parker (1957) and Sweet (1958), the main fea- 
tures of which have been confirmed by many sub- 
sequent calculations (e.g. Biskamp (1993); Uz- 
densky & Kulsrud (2000)). The basic configu- 
ration is sketched in Figure 1, in the coordinate 
system which we will use throughout this paper. 
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Opposing fields ±yBo are pressed together over a 



the inward and outward mass fliixes yields 




Fig. 1. — Configuration for steady state reconnec- 
tion according to Parker-Sweet theory. 

length Ly, leading to a current sheet of width 
centered on 2; = 0. Resistive dissipation of the cur- 
rent annihilates the field. In a steady state, mag- 
netic flux must be brought in at the boundaries at 
the rate at which it is destroyed. The reconnection 
rate is parameterized by the speed at which the 
field is convected inwards, and determining this in- 
flow speed in terms of the geometrical parameter 
Ly and the intrinsic plasma parameters is one of 
the central problems of reconnection theory. 

In Parker-Sweet theory, Uz is determined as fol- 
lows. The electric field E is assumed to be given 
by the standard MHD form of Ohm's law 



u X B J 

E = + -, 

c a 



(1) 



where a is the electrical conductivity. In the 
present problem, E = xE must be constant, be- 
cause the system is in a steady state. We approx- 
imate E by its inductive value — u x B/c every- 
where except in the resistive region, where we ap- 
proximate it by the expression 3 /a. Equating the 
resistive and inductive expressions for E, and us- 
ing Ampere's law to write J as xcBq/{A-kLz), we 
find 



(2) 



where An = c? /{ATia) is the magnetic diffusivity. 

A second relationship between Lz and Uz fol- 
lows from conservation of mass. The inward mass 
flux pUzLy must be balanced by outward mass flux 
in the form of thin jets within the resistive layer. 
Energy conservation arguments suggest that these 
jets travel at the Alfvcn speed ca = Bo/liTrpy^^ , 
implying an outward mass flux pc^Lz- Equating 



Lz — Ly — . 

Ca 



(3) 



Combining equations (2) and (3) yields expres- 
sions for Uz and Lz 
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, AoCA \ c-1/2 
Uz = — — = caS 



Lu 



Lz=(^Y' = L,S-'/\ 

V CA 



where the Lundquist number 5 is defined as 



m 

TA 



Ly CA 



(4) 
(5) 

(6) 



In the interstellar medium, S is generally enor- 
mous (quantitative estimates follow in §2.2), while 
the Alfven crossing time is comparable to the dy- 
namical time. Equation (4) therefore implies that 

reconnection is very slow. 

In this paper, we express the reconnection rate 
in terms of the electric field. In Parker-Sweet re- 
connection, E is given by 



^ ^ UzBq ^ / XnCA\ ^'"^ Bo 
c \ Ly J c 



(7) 



where we have used equation (4). The relations (4) 

and (5) imply that the Ohmic heating rate e^ag 
per unit volume is independent of An 



e-mag — ^2 



47rAnj2 ^ 



47rAn 



-E^ 



(8) 



where the second equality holds only in the re- 
sistive region. The Poynting flux into the region, 
UzBq/{A'k), is proportional to aJj'^^- 

We now consider the problem in weakly ion- 
ized gas. Again, we assume a reversal in yB{z) at 
z = {]. When the charged and neutral matter are 
strongly coupled, regions of magnetic field rever- 
sal can achieve force balance, with neutral pressure 
compensating for the deficit in magnetic pressure. 
Reducing the coupling between neutrals and ions 
leads to ambipolar diffusion, in which case the neu- 
trals lag behind the field lines. Neutral pressure 
support is lost, and the ion pressure - which is 
much smaller than the neutral pressure - is over- 
whelmed by the magnetic pressure squeezing the 
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antiparallel field lines together. The magnetic field 
gradients steepen around the magnetic null, lead- 
ing to high current densities and Ohmic dissipa- 
tion. The basic situation is sketched in Figure 2. 




Fig. 2. — Sketch of magnetic {Pm), ion (Pj) and 
neutral pressure (P„) around the resistive region 
at ^; = 0. The magnetic field strength B would 
increase linearly. The neutral pressure is much 
larger than any other pressure. The magnetic 
pressure squeezes the antiparallel field lines to- 
gether, only resisted by the ion pressure. 

An important simplification of the Parker- 
Sweet problem arises because the flux of ions is not 
conserved. Recombination represents a sink, thus 
eliminating the need for an outflow and making a 
steady state possible in one dimension. Although 
the system could develop 2D or 3D structure - 
and, as we discuss in §6 this may be the only 
way to achieve fast reconnection under certain 
conditions - we treat only the ID case here, be- 
cause it is simplest and most tractable. In the 
spirit of our discussion of Parker- Sweet reconnec- 
tion, we are interested in finding the scaling of 
the steady state inflow speed, or equivalently the 
electric field, with the Ohmic diffusivity An. 

In §2.1, we write down the basic equations 
which govern the system. In §2.2, we nondimen- 
sionalize the problem and define and estimate the 
important basic time scales. In §2.3 we derive the 
reconnection rate assuming ionization equilibrium 
is maintained at all times, and in §2.4 we derive 
scaling relations for the inflow velocity and layer 
thickness when the recombination rate is high but 
flnite. 



2.1. Basic Equations 

BZ95 solved the time dependent layer problem 
for two fluids, one charged and the other neu- 
tral, which are coupled by elastic collisions and 
by ionization and recombination. In the parame- 
ter regime of interest, they found that the problem 
is adequately described by a reduced set of steady 
state equations for the ions. We have generalized 
these equations to allow for a polytropic equation 
of state. Pi oc p] , instead of the isothermal relation 
assumed in BZ95. The dynamics of the neutrals 
are neglected, because the neutrals are not coupled 
to the magnetic field and their pressure profile is 
essentially flat (sec Fig. 3 of BZ95). The ioniza- 
tion fraction is sufficiently small that the neutrals 
can be regarded as an infinite reservoir which can 
be supplemented by recombination with negligi- 
ble effect on its overall properties. The ion flow 
is assumed slow enough that ion inertia can also 
be neglected. As we will see, this imposes a con- 
straint on the width L of the layers that can be 
described by these reduced equations. 

The magnetic induction equation, the ion mo- 
mentum equation, and the ion continuity equation 
are 

|(..P-Aa|i^)=0 (9) 

d 

(PiUi) - RiPnO + RrecPl = (H) 

Table 1 explains the variables. The symbols -R, 

and Rrec denote the ionization rate per unit time 
and the recombination rate per charged particle, 
respectively. Equation (11) is valid for photoion- 
ization or cosmic ray ionization, and for radiative 
or dissociative recombination, but not for surface 
recombination on grains, unless the grains are tied 
to the ions rather than the neutrals. Equation (9) 
can be integrated to yield 

UiB-Xn-^B = -cE, (12) 

and we will always use it in this form. 

Equations (10), (11), and (12) are solved on 
the domain < z < L, with symmetry assumed 
about z = 0. The subscripts "0" denote condi- 
tions at z = L. The boundary values at ^; = are 
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B{0) = 0, Ui{Q) = 0. We also impose the condi- 
tions B{L) = Bq and Pi{L) = pio. We are able 
to impose a fourth condition in a system of three 
first order ordinary differential equations (ODEs) 
by treating E as an eigenvalue. 

For Ao — * the order of equation (12) changes. 
Thus, we expect two physically distinct regions, 
namely a resistive region around z = where 
XnOB/dz dominates, and an inductive region for 
z ^ governed by the inductive term UiB, just 
as occurs in the Parker-Sweet solution. 

2.2. Dimensionless Formulation and Time 

Scales 

As the problem is largely one of time scale ra- 
tios, we rewrite equations (10) - (12) in a dimen- 
sionless form with the relevant time scales as phys- 
ical parameters. We scale the physical variables 
as in Table 1, and introduce the time scales in 
Table 2. We also define the /? parameter 



/3: 



(13) 



The factor of 2 in equation (13) accounts for elec- 
tron pressure, which we assume to be the same as 
the ion pressure. 

We assume ionization equilibrium at the outer 
boundary, such that Tree = Tj. Rewriting the re- 
duced equations as a set of first order ODEs for 
each variable and using equation (13) and Table 2 
we find 



d 



dx 
d_ 
dx 



f3j 
taI 
n r 



1 1 , 
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Tind TA 



OX TA 



rv 



(l-r^) 



V 



d 



r dx 



r. 



(14) 
(15) 
(16) 



Note that the time scales occur only as ratios, 
and that they depend on the length scale L in dif- 
ferent ways: Ti is independent of L, ta is linearly 
proportional to L, and tad and rn are propor- 
tional to LP' . Thus, for example, the ratio t^/tad 
is fixed in terms of the intrinsic parameters of the 
medium, but the ratio of either one to Tj can be 
made larger or smaller by appropriate choice of L. 

Physical realism dictates certain constraints on 
the time scales, which are discussed in the follow- 
ing subsections. 



2.2.1. The resistive time scale tq 

The magnetic difFusivity Xu is related to the 
conductivity perpendicular to the magnetic field 
cr_L by 



An 



47rcr_L 



(17) 



Braginskii (1965) gives (j±_ in terms of the electron 
collision time Te as 



cr_L 



n.e^Te 



(18) 



Electrons collide with both charged and neutral 
species, at rates which we denote by t~^ and t~^ , 
respectively. Combining the effects of collisions of 
both types leads to the composite expression 



TeiTe' 



(19) 



We evaluate Tei from Braginskii (1965) with 
a mean ion charge Z equal to unity, and the 
Coulomb logarithm A set equal to 20. We evaluate 
Ten for H2 as the primary neutral species, using 
the rates given by Draine et al. (1983). The result 
is 



1.4 X 10^2^3/2 



[1 + 1.2 X W-^^T^x-^] 



(20) 



with the ionization fraction Xi = Ue/rin. Electron- 
neutral collisions are represented by the second 
term in square brackets in the denominator of 
equation (20). Such collisions are ignorable in 
most of the environments considered in this pa- 
per. However, under the conditions found in pro- 
toplanetary disks, they dominate over electron-ion 
collisions. 

Equations (17), (18), and (20) combine to give 
the magnetic diffusivity 



2.0 X 10^^ . 



1.2 X lO-^^T^xr^l 



[cm^ s-M 



(21) 



The Ohmic diffusion time tq (see Table 1) is 
then given by 



Tn = 1.54 X 10^^ 



ji3/2 



pc 



[1 + 1.2 X lO-iiT^"^] 



(22) 
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2.2.2. The ambipolar diffusion time scale tad 

The ambipolar diffusivity Xad depends on the 
field strength B, the particle number densities rij, 
n„, the molecular weights /Xj and /z„ (we assume a 
single species of each), the ionization fraction Xi = 
Ui/nn, and the rate coefficient for elastic collisions 
(av). We write 



B 



(23) 



where we have taken {av) = 1.5 x lO'^cm'^s"^ 
(Draine et al. 1983). Thus, the ambipolar diffusion 
time scale is (see Table 1) 



TAD = 9.4 X 10"^ L 



pc 



[yrs]. (24) 



2.2.3. The ionization time scale Ti 



We assume ionization equilibrium at the outer 
boundary so that Tree = Significant departures 
from ionization equilibrium are always in the sense 
of a surplus of ions, so it is the recombination rates 
which matter. We summarize them here. 

Following Draine & Sutin (1987) and McKee 
(1989), we assume that molecular ions are de- 
stroyed primarily by dissociative recombination, 
while atomic ions undergo radiative recombina- 
tion, and also recombine on the surfaces of dust 
grains. Denoting the rate coefficients for these 
three processes by a^r, ccgr, and a^r, respectively, 
we write the total recombination rate Rrec as 

Rrec = ifiirriHy^ [(t>adr + (1 - (p)arr + CtgrX~^] , 

(25) 

where ^ = for atomic ions and (j) = 1 for molec- 
ular ions. The recombination time including all 
three processes is then given by 



'^rec (-^recPi) 



(26) 



Equation (25) is valid when the grain density 
follows the neutral density. As such, it is incon- 
sistent with equation (11), because the rate of re- 
combination on negatively charged grains is then 
linear in the ion density rather than quadratic. 



We have not included this possibility in our calcu- 
lations, but we comment on it in §6. It is possible 
that at least some of the grains follow the ions; 
sec Ciolek & Mouschovias (1993) for a discussion 
of grain-gas coupling. 

Taking expressions for the a coefficients from 
Draine & Sutin (1987), we have 

7.2 X 10-- ^^^^ 



O^dr 



rpl/2 

1.4 X iO"-^r--''/4 

2.9 X 10-12 / 3A 



gr 



3/2 



(28) 
(29) 



all in units of cm^s-^. In equations (29), amin is 
the minimum grain size, with the size spectrum 
assumed to follow the MRN power law (Mathis et 
al. 1977). For the numbers given in Table 3, we 
used amin = 3 A. Although this value of amin is 
uncertain, our final results are rather insensitive 
to it, because of the relatively large coefficients of 
radiative and dissociative recombination. 

2.2.4. The Alfven crossing time ta 

The non-dimensionalization leaves us with a 

"crossing time" in the problem. There are vari- 
ous ways to define this crossing time. As it is not 
clear what the velocity at the outer boundary will 
be (see below as well) , but values for the magnetic 
field Bq and the density piQ are given, we scale the 
velocity in terms of the ion Alfven speed cai at 
the outer boundary .t = 1, thus introducing the 
Alfven crossing time of Table 1 . We have 



CAi = 2.2 X 10 
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B 



The corresponding crossing time is 



cm s 



(30) 



TA = 0.46 Lpc- 



B 



[yrs]. (31) 



2.2.5. The inductive time scale Tind 

The inductive time Ti„d plays a special role. 
Unlike the other time scales, it is not a free pa- 
rameter, but acts as an eigenvalue of the set of 
ODEs (14 - 16), allowing us to impose four in- 
stead of three boundary conditions. We define the 
inductive time by the scaled electric field 

(32) 
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2.2.6. The length scale L 

Wc have not specified the length scale L in 
equations (22) and (24). It is controlled in part 
by the global setup of the problem, but is gov- 
erned by certain constraints. We require L to be 
small enough that the neutrals are decoupled from 
the magnetic field, but large enough that the ion- 
neutral collision time is shorter than the local 
ion dynamical time. In the subsequent discussion, 
we leave L unspecified, in order to keep the argu- 
ment as general as possible. In §5 we discuss the 
physical systems where we believe the combina- 
tion of ambipolar diffusion and reconnection to be 
relevant. Table 3 contains characteristic numbers 
for diffuse clouds, dense clouds, cores and proto- 
planetary disks. The numbers demonstrate the 
disparity of time scales in the problem. 

2.3. Zero Ion Pressure Gradient 

If there is no ion pressure gradient, the dynam- 
ics are extremely simple, and serve as a reference 
to which more general solutions can be compared. 
This limit is achieved in the case /? = and/or in 
the case = 0. We refer to this solution as the 
pure ambipolar diffusion (AD) solution. 

According to equation (15) 



V = b — b = VAD- 

TAD OX 



(33) 



With equation (33), the magnetic induction equa- 
tion (14) can be integrated to yield 



1 



.53 + _5 = 

TQ, 



c. 



(34) 



Since & = at a; = 0, the constant of integration 
C must be zero. Using the boundary condition 
6(1) = 1 and assuming m tad, we get for the 
scaled electric field its AD-solution, namely 



cE 



1 



^TAD 



(35) 



Thus, under the conditions of pure AD, with 
^n/^AD ^ Ij -E' is independent of Xq, so recon- 
nection is fast. 

Equation (34) predicts that the solution has a 
boundary layer near x = 0. Far from the origin. 



JAD 



VAD 



,1/3 



^TAD 



(36) 



Near the origin, 

bAD ' 
VAD '■ 



-X 



ta 

TAD 



TQ 



3ta 



D 



Equation (36) holds for x > xad, where 



XAD = 



D 



3ta 
Tn 



3/2 



(37) 



(38) 



For X < XAD, the solution is given by equa- 
tion (37). 

As equations (36) and (37) show, the flow speed 
increases toward the origin, then plunges to zero. 
Thus, in the linear region, the ions are subject to 
tremendous compression, raising the possibility of 
strong deceleration of the flow by the ion pressure 
gradient. In order to see how this limits the rate of 
reconnection, we consider the role of ion pressure. 

2.4. The Effect of Ion Pressure 

As long as <C 1, and tq <^ tad, ion pressure 
and resistivity are important only within a bound- 
ary layer defined by a; < a;;. Regularity requires 
V oc X as X ^ 0, and so we introduce a dynamical 
time Td such that 



TA 

V = X, 

Td 



(39) 



for X < xi. 

Equation (39) permits us to solve the continuity 
equation (16) for the central density r(0) = r™. 
The result is 



2rd 



1+ - 

\2Td 



-1 1/2 



(40) 



We expect > 1, for large tq even 1, so 

that we can write 



Td 



(41) 



We now take the opportunity to correct an er- 
ror in BZ95. Symmetry considerations require 
that the Taylor series for r{x) contains only even 
powers of x: r{x) = + r2x'^ + while the 
Taylor series for v{x) contains only even powers: 
v{x) = —TAx/Td — v^x^ — .... BZ95 solved for r2 
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without including W3. This is incorrect, although 
their conclusions as to the reconnection rate are 
not affected in the parameter range they consid- 
ered. 

We derived equations (36) and (37) by setting 
the constant of integration C in equation (34) 
equal to zero. Now, however, we are not following 
this solution all the way to the origin. This allows 
us to introduce a parameter a;* such that 



b = 



X + .X, 
1 + 



1/3 



V = —- 



The scaled electric field is then 



£ = 



3ta£)(1 +a;*)' 



(42) 



(43) 



so X* is effectively a parameterization of £. For 
a;* = we recover the ambipolar diffusion solution 
equation (35), and this solution is approximately 
valid as long as < 1. 

Our description of the system depends on three 
parameters which are so far unknown. The first, 
Xi , is the location of the transition between the in- 
ner, resistivity and pressure dominated boundary 
layer and the outer, inductive region. The sec- 
ond, Td, is the dynamical time in the inner layer. 
The third, , measures the departure of the outer 
solution from the pure AD solution. We can de- 
termine these three parameters, thereby patching 
together a global solution, by applying three con- 
ditions which must connect the inner and outer 
domains. 

First, we assume that at x = the solution 
changes from inductively dominated to resistively 

dominated. Equating the left hand side of equa- 
tion (14) to the second term on the right hand 
side, and using equation (39), we find 



Xl= [ — 



1/2 



(44) 



As the second condition, we assume continuity 
of total (plasma plus magnetic) pressure at xi , and 
also that magnetic pressure dominates for x > Xi 
while plasma pressure dominates for x < x;. Our 
picture is that the plasma pressure profile inside 
X = X; is quite flat. These assumptions together, 



along with equations (15) and (42), lead to 



y ' m 



XI + X* 

1 -I- X* 



2/3 



(45) 



Finally, we postulate mass flux conservation 
across xi. Using equations (39) and (42), we have 



Xl 



^rf 3taz3(x;+X*)^/^(1 + X*)^/^" 



(46) 



Eliminating in equation (41) with equa- 
tion (44) yields 



Taxf 



(47) 



Equation (47) together with equations (46) and 
(45) leads to a pair of equations for xi and x*: 

x^i\xi+x.f/\l + x.)--"^ = /3(^)'(48) 



xf(x(+x.)-i/3(l+x.)-2/3 = 5^ 



(49) 



It turns out that x^/x* <^ 1 for all solutions 

of equations (48) and (49). Therefore, we can ap- 
proximate equations (48) and (49) by 

xf^x^/'(l + x*)-2/3 = (50) 



(51) 



Raising equation (50) to the power 3/(27) and di- 
viding by equation (51) yields an equation for x* 
alone 



3±1 2-,-3 , ,„ 

x/^ (l + x.)~ =^'/', 



where 



Z = /33/' 



TjTQ 



(52) 



(53) 



According to this analysis, Z is the fundamen- 
tal parameter which controls the properties of the 
layer. When Z <C 1, x* <C 1 as well, and equa- 
tion (52) has the approximate solution 



X* W < 1. 



(54) 
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When Z ^ 1, Xi, ^ 1, and equation (52) yields 
w Z^/'^ > 1. (55) 

Equations (54) and (55) agree closely when 
7 = 5/3, in which case the exponent in equa- 
tion (54) is 15/28 « 0.54. However, it is useful for 
comparison with numerical work to find solutions 
of equations (50) and (51) which are valid for all 
a;* ^ X;. Dividing equation (50) by equation (51), 
we get 



27-3 
Xt x^ 



fill 

3 



7-1 2--7 
TAD 



Q. 



We can express a;, as 

and with equation (50) arrive at 



l + Qx 



3-27 



/3 



3/2 



(56) 
(57) 

(58) 



Equation (58) can be solved numerically for xi, 
with X* determined from equation (57). The re- 
sults arc shown in Figure 3 for parameters taken 
from the series of numerical models A listed in 
Table 4. 

It is also useful to have approximate analytical 
expressions for the boundary layer width X; in the 
limiting cases Z ^1, Z ^ 1. With equations (50) 
and (54) we find 



^7+2^2 



1/(6+27) 



XI = 9/3- 



' AD 



^7+4 



(59) 



for Z <C 1. Using equations (50) and (55) we have 

1/2 



Tn 



(60) 



for Z > 1. 

Finally, wc use equation (43) to derive the scal- 
ing of the electric field E with magnetic diffusivity 
An- For simplicity we assume x* = Z^/^ for all Z. 
We find 



dlnS 



Z 



dlnXn 2{1 + VZ)' 



(61) 



For Z ^ 1, £ (X A^/^, as in Parker-Sweet recon- 
nection, and for Z <C 1 we expect £ oc X^, i.e. 
"fast" reconnection. 



4.0 
m 0.0 



c 
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0.52 



0.4f 



-4.0 I*** 
0.56 



o .A34 
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□ .^32 



10 12 14 K 



log T5, 

Fig. 3. — X* against tq for the general solution 
of eqn. (58). The exponent change connects the 
regimes x; <C x* <C 1 and x; <C 1 <C x*. 

3. Numerical Method 

We seek a numerical solution of the three cou- 
pled ODEs (14) - (16). As we mentioned below 
equation (12), we require that the three depen- 
dent variables &, r and v satisfy four boundary 
conditions (BCs): 



6(0) = 
u(0) = 
6(1) = 1 
r(l) = l 



(62) 
(63) 
(64) 
(65) 



The first two conditions follow from symmetry ar- 
guments, and the second two from our character- 
ization of the equilibrium state far from the re- 
sistive layer. Note that we do not specify w(l), 
leaving it free to adjust to the varying resistivity. 

Equations (14) - (16) together with condi- 
tions (62) - (65) constitute a two-point boundary 
value problem which is also an eigenvalue prob- 
lem for £. We implement the constraint on £ by 
introducing a fourth equation 



—£ - ^ 

dx dx Tind 



0. 



(66) 



We solve the problem using a combination of 
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shooting and relaxation techniques (Press et al. 
1992) 

The shooting method is based on treating the 
problem as an initial-value problem with bound- 
ary conditions implemented only at a; = 0. Of A'^ 
boundary conditions, let rii be imposed at a; = 
and 712 = N — m BCs at x = 1. Thus, for the 
initial- value problem, we have n2 freely specifiable 
parameters at a; = 0. Outward integration results 
in values of r, v and 6 at a; = 1 which probably do 
not fit the corresponding BCs. We zero the differ- 
ences via Newton-Raphson root finding and thus 
find corrections to the n2 free values at a: = 0. The 
process is repeated until the differences at a; = 1 
are below a given threshold. 

The relaxation method uses a finite-difference 
approximation of the ODEs on a grid with given 
number of grid points. It too is based on a New- 
ton root-finding scheme (Press et al. 1992). We 
use the scheme in its simplest form, coupling two 
neighboring grid points. Higher order coupling is 
possible, but renders the method much more com- 
plex. 

The relaxation method requires an initial guess 
over the full domain for all variables. Thus, we 
combined both solvers, using the shooting method 
to obtain an initial guess for the relaxation. 

The length scales and time scales in this prob- 
lem vary over large ranges. For example, the ratio 
Ta/TAD ~ 10^° for physically realistic parameters. 
Therefore, our set of ODEs is very stiff. Neither 
the standard shooting method nor the relaxation 
method can solve the equations for such param- 
eters. We resolve the problem with a threefold 
approach. 

First - as the most obvious measure - we use 
a non-uniform grid. The standard solution would 
be a logarithmic grid, however, the normalization 
x{N) = 1 comes naturally with an exponential 
grid 

x{k) = e^^'^-^K (67) 

N is the number of grid points, k the grid point 
index, and 6 a compression factor. The boundary 
conditions (eq. 62-65) require 6(0) = and v{0) = 
0. As wc cannot reach x = with equation (67), 
we interpolate linearly between a cutoff Xcut and 
a; = 0. We choose Xcut such that we do not lose 
any resolution in physically interesting regions. 



Second, we approach physically realistic param- 
eters from an initial guess by taking small steps 
in parameter space. For the A series of runs, we 
began with unrcalistically small values of tq and 
advanced slowly to more realistic values. 

Finally, we found that at large tq or small (3 the 
density and velocity gradients are so steep that we 
can resolve them with the exponential grid only 
by using an enormous number of grid points, far 
more than are needed over most of the domain. 
Accordingly, we introduced a simple version of 
adaptive mesh refinement. We refine on the ve- 
locity difference between two adjacent grid points, 
Aiifc = \vk-Vk-i\- If Aufc > Av+ for ki < k < k2, 
the number of grid points in this region is multi- 
plied by n+ = 8. If Avk < Aw_ for kx < k < k^ 
under the condition k2 — ki > ri-, the number of 
grid points in this region is divided by n_ = 8. 
We do not allow refinement or coarsening at the 
lower and upper boundaries, which in any case is 
unnecessary. Moreover, refinement is not allowed 
for 

— < Emacft, (68) 

Xk-1 

where emach is the machine accuracy. The a;-grid is 
interpolated linearly on the new grid points. The 
number of gridpoints N typically varies between 
lO'^ and 3 x 10^ 

Convergence to a stable solution depends cru- 
cially on the overall error criterion. We enforce 
the fractional error econv = Ay/y < 10""* for all 
gridpoints, where y is any of the physical variables 
b, r, V and £. In order to ensure that we find the 
correct solution for each parameter adaption itera- 
tion, the successive change of a parameter, e.g. tq, 
must lead to differences in the results larger than 
the convergence error ^conv Thus, econv needs to 
be adapted with decreasing parameter adaption 
step size. 

We have confirmed that the solution remains 
independent of N, as long as N is large enough 
for the relaxation method to converge. The criti- 
cal region is near the edge of the boundary layer. 
Solutions for different A^ are identical within the 
convergence error econv The method develops 
oscillations in the density aromid x < 1 with 
r/r exact < ^conv, but the overall solution is not 
affected. 
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4. Results 

We first describe the basic properties of the so- 
lutions and discuss the differences between simula- 
tions and analytical predictions. Covering a range 
of 8 orders of magnitude in tq, we then can in- 
fer scaling laws from the simulations and compare 
them to the predictions. 

In order to investigate the dependence of the 
scaling relations on and /?, wc ran models A 
within four parameter sets, which are listed in Ta- 
ble 4. 

4.1. Basic Properties of the Solutions 

Figure 4 shows the full domain for all mod- 
els of type ^ at TQ = 10^^. Gas is moving in- 
wards, i.e. to small x, at approximately constant 
velocity, until it is brought nearly to a halt by 
the rapidly increasing plasma pressure at xi. The 
magnetic field and the velocity are approximately 
linear within the boundary layer, as expected from 
equations (37), but essentially flat in the ambipo- 
lar or inductive regime. From equations (36) we 
would expect b cx x^^^ and v oc x~^^^. Only 
for model A42 do we notice a slight hint of such 
behavior. However, in Figure 5, the overplotted 
outer solutions of equations (36) fit perfectly to 
the numerical results. The fractional differences 
between the numerical and analytical solutions at 
the transition point between the resistive and in- 
ductive region amounts to 



bAD — bsim 




VAD — Vsim 


bAD 




VAD 



10" 



(69) 



where the subscripts AD and sim denote the outer 
(eqs. 36) and simulation's solution. We realize 
that the solutions are so flat because, aside from 
a small region around x = \, we have a;* ^ a; for 
all models except A42 (see Figure 3). 

The analytical prediction of the reconnection 
rate (§2.4) is derived assuming pressure balance 
between magnetic and ion pressure (eq. 45) and 
conservation of mass flux across the transition 
layer (eq. 46). Figure 6 confirms these assump- 
tions. Note that the mass flux scale in the lower 
panel is linear, while the velocity scale in Fig- 
ure 4 is logarithmic. While the velocity jumps over 
nearly two orders of magnitude within Aa; « 10~^, 





m -2 
o 

-4 
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-4 
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log X 

Fig. 4. — Magnetic field 6, density r and velocity 
u at rn = 10^^. The resistive and ambipolar diffu- 
sion regimes are separated by a sharp transition. 

the mass flux changes at most by a factor of nearly 
4 over at least one order of magnitude in a;. 

The analysis in §2.4 takes the inductive and 
resistive contributions Emd and f^es to £ to be 
equal at xi, and also uses the assumed linearity 
of V and h with x. Figure 7 shows xi vs tq and 
d{lnxi)/dlnTn vs tq for the A simulations and 
for the corresponding analytical predictions. The 
values of xi in the simulations can be measured 
accurately and unambiguously to three significant 
digits from the first or second derivatives of any 
of the variables b, r and v. The predicted val- 
ues of xi are derived from equation (58). The 
upper panel of Figure 7 shows a persistent offset 
between the simulated and predicted values a;f"" 
and xf^'^, in the sense that x^*"* is always less than 
xf^'^ , sometimes by more than an order of magni- 
tude. As the lower panel of Figure 7 makes clear, 
the slopes of the xi{tq) relations differ as well, al- 
though d{liixi)/d{\nTa) computed from the simu- 
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Fig. 5. — Magnetic field b and velocity t; at rn = 
10^^. Thick lines: numerical solution, thin lines: 
outer solution according to equations (36). 

lations converges to the predicted value of 0.5 for 
large tq. 

In order to understand this discrepancy, we 
note that as long as v is given by equation (39), 



2 ^ ^ind Td d\nb 



(70) 



where all quantities on the right hand side are 
evaluated at xi. Equation (70) reduces to equa- 
tion (44) for £ind/£res = 1 and rf(ln6)/rf(lna;) = 1. 

Examination of d{h\b) / d{\n. x) in our solutions 
shows that it drops sharply as a; — > xi. With 
the diminished importance of the resistive term, 
Sind/^res = 1 at actually exceeds unity for most 
runs. However, db/dx is steeper in the inner parts 
of the layer than it would be if d(ln b) / d{h\ x) were 
equal to one throughout. This cfTcctivcly makes 
the layer thinner than equation (44) predicts. 

4.2. Reconnection Rates 

The electric field determines the magnetic re- 
connection rate (eq. 7), and thus the magnetic 
energy dissipation rate per volume e-mag (eq- 8). 
BZ95 implicitly argued for E oc A^^, or e-mag oc 
Aq^, meaning fast reconnection. They used 40 < 
to/t^d < 4 X 10^, resulting in a range 4 x 10~^ < 
Z < A X 10^. Thus, they were well within the 
fast reconnection regime. However, they assumed 
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Fig. 6. — Upper panel: Magnetic and ion pres- 
sure, and their sum for all models at to = 10^^. 
Magnetic and ion pressure add up to a flat curve 
across xi. For x > 10~^, the magnetic pres- 
sure is balanced by the third term in equation 10. 
Lower panel: Mass flux for all models as above. 
Compared to the corresponding change in veloc- 
ity (up to two orders of magnitude), the mass flux 
is nearly constant at x;. 

isothermality, leading to higher ion densities in the 
resistive layer. 

In the previous section we saw that the scal- 
ing for xi changes at a;* « 1. As Figure 8 
demonstrates, this extends to the electric field 
scaling. The 8 orders of magnitude in tq cover 
the "fast" regime with £ independent of tq and 



the "slow" regime with S (x 



-1/2 



1 /2 

£ oc Ao (lower panel of Fig 



or equivalently 
i^iuwci pciiici ui 1. ig,. oy. Both regimes 
exist. Whether they are both equally realistic, we 
will discuss in §5. 

The lower panel of Figure 8 shows that the 

scaling of £ with An approaches the predicted 

1/2 

£ oc Aq at large rn. However, there are two dis- 
crepancies between the simulations and the theory. 
First, the value of £ predicted by the analytical 
theory is always larger than the value found in the 
simulations. Second, the lower panel shows that 
the predicted value of to at which the reconnec- 
tion rate changes from fast to slow is larger than 
the simulated value. 
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Fig. 7. — Upper panel: Resistive scale xi against 
resistive time tq, predicted (lines) according to 
equations 50 and 51, and determined from the sim- 
ulations (symbols). Lower panel: Exponent q of 
xi oc Tq. The dotted line denotes q = —0.5, the 
theoretical value for a;* 1. 



In the previous subsection, we asserted that two 
of the key assumptions of the analytical theory - 
overall pressure balance across the layer, and con- 
servation of mass across the transition - are vali- 
dated by the simulations. The third assumption, 
equation (44), is not as well supported because of 
the flattening of b with x in the outer parts of the 
layer. This leads to an overestimate of xi relative 
to the simulations. 

Much of the discrepancy between £sim and 
can be traced to the overestimate of xi in the an- 
alytical theory (fig. 9). Symbols denote £smi/£pre 
for the A sequence models. Lines stand for the 
same quantity with the predicted xi replaced by 
Xi from the simulations. It is clear that the agree- 
ment between theory and simulation is much bet- 
ter in the lower panel. 

Theory and simulation agree significantly bet- 
ter for series A2>2 and ^42 than for series ^34 
and >144 (figs. 8 and 9). In the first two series of 
models, the inequality 



n «; TAD < Tn 



(71) 



O 




is well satisfied, while for the other two series it is 
not. 



8 10 12 

log 

Fig. 8. — Electric field £ against resistive time rn 
for all models. Lines represent theoretical predic- 
tions, symbols stand for simulations. The lower 
panel shows the exponent p as in f a Tq. 

In summary, we see a reasonable degree of con- 
sistency between the simulations and the analyti- 
cal theory. Based on their agreement, we are confi- 
dent that the theory identifies the important phys- 
ical effects, while the simulations provide numeri- 
cal factors which correct the theory. 

According to the theory, the reconnection rate 
is determined by the parameter Z defined in equa- 
tion (53) . Figure 10 shows £ and dln£/d\n.Z plot- 
ted against Z for all the A series simulations (the 
usual symbols) and for the theory, shown as a solid 
line. The first observation is that Z is indeed a 
good parameter for ordering the results of the sim- 
ulations. While the simulations do not fall onto a 
single curve because of the effects discussed above 
(e.g. eq. 71), they lie quite close together. The 
analytical predictions for different models are of 
course indiscernible. We see also that reconnection 
slows down at a smaller value of Z than predicted 
by the theory. Empirically, we suggest Z w 0.01 
rather than Z = 1 as the dividing line. We iden- 
tified the discrepancies between predictions and 
simulations as being caused by the over-idealized 
linearity of b and v for x < xi. Thus, we will use 
the analytical predictions for the subsequent dis- 
cussion, but correct them by replacing Z by 100 Z. 
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Fig. 9. — Ratio of simulated over predicted electric 
field £ against resistive time tq for all models. The 
agreement between theory and simulation deterio- 
rates for large Symbols denote £sim/£pre, lines 
stand for the same quantity with the predicted xi 
replaced by xi from the simulations. 

5. Physical Interpretation 

We now apply our rc;siilts to weakly ionized as- 
trophysical environments. We begin by discussing 
the choice of length scale L. Then we derive the 
rcconncction rate as a function of L for some astro- 
physical systems of interest. We estimate the rate 
of Ohmic heating and show that it is extremely 
high within the resistive layers. Finally, we discuss 
the implications for the ISM and for protoplane- 
tary disks. 

5.1. The Choice of L 

One crucial parameter has yet to be defined: 
the length scale L of the domain (see §2.2.6). The 
scaling parameter Z - and thus the rcconncction 
rate - depends on L as L~'^. Thus, we predict 
faster reconnection in smaller systems, which is 
intuitively plausible. 

In order for the set of equations (9) - (11) to 
be valid, the ions and neutrals must move inde- 
pendently of one another. To estimate the length 
scale at which this occurs, we consider small per- 
turbations of the two fluid system having wave 
number k perpendicular to the magnetic field. In 
the limit fc ^ 0, the ions and neutrals arc almost 
perfectly coupled. As k increases, the coupling de- 
creases, until I {vi — Vn)/vi I 1 at a wavenumber 
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Fig. 10. — Electric field £ against scaling parame- 
ter Z as defined in equation (53). 



ko which is given by 



CD 



2Vn 
C-An 



1 + 



1/2 



(72) 



Thus, we require L < 27r /kn = Lad- For small 
Csn/cArn cquation (72) reverts to the usual re- 
sult for hydromagnetic waves (e.g Ferriere et al. 
(1988)). If Csn/cAn is large, as we might expect in 
the vicinity of a neutral sheet, the critical scale is 
smaller by a factor of CAn/csn- 

Lad can also be defined in terms of the am- 
bipolar Reynolds number Rad (Zweibel & Bran- 
denburg 1997; Klessen et al. 2000) 



Rad = 



Lv 

^AD ' 



(73) 



where Xad is defined in equation (23) and v is the 
speed of the neutrals. Then, Lad is the value of 
L for which Rad = 1. If we replace v by the mag- 
netosonic speed {c\^ + c^„)^/^, then Lad ~ kj)^- 
However, equation (73) also allows for the possi- 
bility that the system is dominated by supcrmag- 
netosonic fiows. It is possible that L could be less 
than Lad if it were defined by a geometrical prop- 
erty of the system, such as the thickness of an ac- 
cretion disk. However, L cannot be greater than 
Lad if the equations we solved are to be valid. 
The requirement that the system be collisional 



14 



gives a lower limit for L, namely 

L Lcol = Csi/l^iri') • 



(74) 



where = IJ'nnn{(^'o) / Hn) is the ion-neutral 
colhsion frequency. 

5.2. Reconnection Rates in Physical Envi- 
ronments 

Table 3 lists representative values of physi- 
cal parameters in four environments: (a) diffuse 
clouds, (b) dense molecular clouds, (c) molecular 
cores and (d) protoplanetary disks. In Figure 11 
we plot the incoming ion flow speed - the recon- 
nection rate - in units of the neutral Alfven speed 
CAn against the scale length. 
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Fig. 11. — Scaling parameter a;, and reconnection 
rate Ui in units of the neutral Alfven speed cau 
against domain scale for the parameter sets (a), 
(b), (c) and (d) as defined in Table 3. The ver- 
tical lines denote the constraints on L discussed 
in §5.1, and right (left) pointing arrows stand 
for lower (upper) thresholds, according to equa- 
tions (73) and (74). We use the empirical formula 
= 10^1/2. 

At small L. the ciirves reach a saturated, max- 
imum value Uis' 

Using eqns. (42) and the upper 
panel of Figure 11, we see that this saturated ve- 
locity is 



An 



1/2 



(75) 



At large L, eqns (42) show that Ui is given approx- 
imately by Xad/'^L. However, only the part of the 
domain with Lcol < L < Lad is accessible for any 
given environment. These limits are marked on 
each of the curves. 

Figure 11 shows that magnetic fields can merge 
quite rapidly under realistic astrophysical condi- 
tions. Note that the reconnection speed Uis, which 
is the maximum reconnection speed as a function 
of L with all other parameters held fixed, actu- 
ally corresponds to slow reconnection in the sense 

1/2 

we originally defined it: Uis (x Xq . However, Uis 
is generally much faster than the Parker-Sweet re- 
connection speed Uz defined in equation (4). Using 
equations (6) and (75), we see that 



TA 



(76) 



where ta is to be computed using the total mass 
density. The relative rapidity of magnetic merging 
comes about because the recombination time is 
short compared to the flow time along the neutral 
sheet, and the plasma (3 is very small. 

5.3. Ohmic Heating 

If any aspect of magnetic merging of the type 

described here is observable, it is likely to be re- 
lated to Ohmic heating. In a steady state, the 
flux of magnetic energy into the layer is balanced 
by radiation, resulting in a radiative flux of 



Fr = 2u., 



Bl 



(77) 



The spectrum of the emergent radiation de- 
pends on the temperature in the layer. While a 
full thermal equilibrium model is beyond the scope 
of this paper, it is useful to calculate the Ohmic 
heating rate. At the midplane, according to equa- 
tions (8) and (43), the energy dissipation rate per 
unit volume is 



^maa — 



AlTTifi^h (l+x*)2 



(78) 



Recalling that x* ~ Z^/^, we see that equa- 
tion (78) predicts that Cmag oc L"^ for a;* 4; 1 
and reaches a saturated value of BQ/{A'KTi (3^/"') 
for ;» 1. We have evaluated e^ag for the mod- 
els in series A. It is nearly constant at the mid- 
plane value throughout the inner boundary layer. 
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and plunges to very small values outside it. Fig- 
ure 12 shows the boundary layer magnitudes of 
Cmog as a function of L for the four environments 
listed in Table 3, with the permitted range of L 
values indicated in each case. The dependence on 
L follows the prediction of equation (78). 
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Fig. 12. — Ohmic heating rate Cmag against scale 
length L for the parameter sets (a) . (b) , (c) and (d) 
as defined in Table 3. Physically relevant regimes 
are denoted by vertical lines as in Figure 11. We 
used the empirical result a;* = lOZ^/^. 

The saturated values of e,„ag in all these en- 
vironments are extremely large in comparison to 
the radiative cooling rates. For example, in dif- 
fuse clouds [case (a)], the maximum cooling rate 
by the 158/im fine structure line of C II is 8 x 
10~^^ ergs cm~^ s~^. This implies that the gas is 
heated to very high temperatures, and perhaps is 
even thermally ionized. This must slow down the 
reconnection rate at smaller values of tq than pre- 
dicted by our theory. We could very roughly ac- 
count for the increase in plasma pressure by choos- 
ing 7 > 5/3. The theory predicts that as the equa- 
tion of state stiffens, the reconnection rate slows 
down and the boundary layer thickens. 

The high current density in the boundary layer 
might destabilize high frequency waves which 
would provide anomalous resistivity through 
wave-particle interactions and increase the merg- 
ing rate. Following the early suggestion by Krall 
& Liewer (1971), attention focussed on the lower 
hybrid drift instability as the mechanism for gener- 
ating the waves. Although it now appears unlikely 
that drift waves themselves are a source of anoma- 
lous resistivity, but they may play a role in the 



growth of Kelvin-Helmholtz instabilities, which 
change the topology of the reconnection region 
(Lapenta & Brackbill 2002). Additional sources 
of small scale turbulence in null layers could also 
contribute to anomalous resistivity (Ozaki et. al. 
1996). 

Although the heating is intense, the narrowness 
of the layers is a formidable barrier to their di- 
rect detectability. The analytical formulation pre- 
dicts that the dimensional layer width, Lxi, scales 
as L^/^^ for Z <C 1 and is independent of L for 
Z > 1; Lxi ~ (/Ji/'^r, Aa)i/2. The layer widths 
are plotted in Figure 13 for the four environments 
listed in Table 3. 
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Fig. 13. — Resistive layer width Lxi against scale 
length L for the parameter sets (a), (b), (c) and 
(d) as defined in Table 3. At « 1 (see fig. 11), 
the dependence of xi on L changes from flat to 
mcreasing. We took X* = 10ZV2 



5.4. Reconnection in the ISM and Proto- 
planetary Disks 

Figure 11 allows us to read off the range of 
Icngthscales over which our theory applies, and 
rates of reconnection predicted by the theory, in a 
given environment. For the diffuse clouds (a), the 
reconnection rate saturates at a sub-Alfvenic value 
of Uis- Within the permissible range of length 
scales - denoted by the vertical solid lines and 
inward-pointing arrows -, Ui has reached its satu- 
rated value Uis- Reconnection can occur only if the 
system survives for a few tens of Alfven crossing 
times. 

The situation is different for the environments 
(b) and (c), namely molecular clouds and cloud 
cores. In both cases, the saturated reconnection 
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rate Uis is super- Alfvenic by several orders of mag- 
nitude, and the permissible scale length range ex- 
tends well into the saturated regime. Even in 
the unsaturated regime, it is possible to merge 
at roughly CAn if i is a few hundredths to a few 
thousandths of a parsec. Such scales can be im- 
aged directly, and arc known to support density 
structure. Therefore, reconnection is a viable pro- 
cess for the annihilation of opposing flux systems 
in molecular clouds. 

Protoplanetary disks are much smaller than in- 
terstellar clouds, so curve (d) of Figure for L < 
10~^ pc. The physically permissible scale lengths 
range from Ri lO"'^ pc down to a few meters. 
The inflow velocity saturates only at the small- 
est scales, but the flow velocity is super-Alfvenic 
over approximately 7 orders of magnitude in scale 
length L. The heating provided by fast reconnec- 
tion (Fig. 12) points into the direction of models 
for chondrule formation as discussed by e.g. Levy 
& Araki (1989), Cameron (1995) and Joung et al. 
(2002). 

6. Summary and Discussion 

In this paper, we have determined the rate at 
which magnetic flux can be supplied to a weakly 
ionized, neutral layer where it is resistivcly de- 
stroyed under steady state conditions. Through 
an extensive parameter study, covering 8 orders of 
magnitude in resistivity, we were able to find the 
regimes of slow and fast reconnection, and iden- 
tify the parameter Z (eq. 53) which distinguishes 
between them. Our simple analytical model and 
the simulations agree qualitatively. Quantitative 
agreement follows directly from accounting for 
the over-idealized linearity of the analytical solu- 
tion. Thus, we are confident that we can apply 
our model to astrophysical environments, provided 
that no additional effects intervene. 

Our results generalize and extend BZ95. In 
§2.4, we pointed out an error in the formulation 
in that paper. The analysis in BZ95 predicts re- 
connection at a rate independent of Ao as long as 
Z < p, while the analysis here predicts that the 
rate is independent of Ao as long as Z < 1. 

Our model is the simplest one possible. The 
flow is one dimensional, and the magnetic field 
vanishes at the midplane. Thus, it is worth ask- 
ing whether the model can be generalized. Two 



dimensional flow, as in the Parker-Sweet model 
of reconnection, provides another mechanism for 
removing ions from the resistive layer. However, 
we already know that this leads to a slow re- 
connection rate. Fast reconnection in our model 
is achieved through a short recombination time, 
which removes ions more rapidly than outflow at 
the Alfven speed. The 2D problem was consid- 
ered by Vishniac & Lazarian (1999). Their result 
for the reconnection rate, if Ohmic heating is neg- 
ligible, is slower than ours by a factor (3~^. 

If we included a third component of magnetic 
fleld, which did not vanish at the neutral line, the 
modifications would be more substantial. Mag- 
netic pressure from this third component would 
reduce the flow speed and prevent steepening of 
the profile. In the case of a dominant third com- 
ponent, we should recover the theory of tearing 
modes in weakly ionized gas (Zweibel 1989), in 
which the reconnection rate is enhanced only by a 
factor of {p/piY^^, and thus remains slow. 

Reconnection is also somewhat slowed down 
in situations where recombination proceeds on 
grains, and the grains are tied to the neutrals. 
An analysis similar to that in §2.4 shows that 
X, « /3-i/27^i/2 for ^ 1, 

Bearing these caveats in mind, fast reconnec- 
tion - i.e. a reconnection rate independent of An 
- is possible within the model presented. The con- 
ditions which must be met are (i) a short ioniza- 
tion time, guaranteeing the fast removal of ions 
in the resistive layer and thus preventing a pres- 
sure build-up, and (ii) a negligible ion pressure in 
comparison to neutral and magnetic pressure. 

As we saw in Figure (11), the maximum merg- 
ing speed is actually "slow" in the sense that it 

1 /2 

depends on An as . Quantitatively, however, 
it is quite fast: super-Alfvenic for the molecular 
cloud and protoplanetary disk environments which 
we considered (cases (b)-(d) of Table 3), and sub- 
Alfvenic by less than two orders of magnitude for 
diffuse clouds [case(a)]. We showed that at length- 
scales less than the ion-neutral decoupling length 
Lad defined by equation (73), reconnection is in- 
deed in the super-Alfvenic regime. In protoplane- 
tary disks (case (d)), fast reconnection might be a 
good candidate to supply sufficient heat for chon- 
drule formation. Fundamentally, reconnection is 
rapid as long as the plasma (3 is low and the re- 
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combination time of ions is much less than the flow 
time at the Alfven speed. 

An extension of the model to higher dimensions 
and including time dependence seems desirable, in 
order to test whether the process is indeed influ- 
ential under more realistic conditions. 
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Table 1 

Physical and dimensionless variables 



quantity 


physical 


dimensionless 


length 


z 


X = z/ L 


ion density 


Pi 


r = Pi/Pio 


ion velocity 


Ui 


V = Ui/cAi 


magnetic field 


B 


b = B/Bo 


electric field 


E 


£ = cE/{LBo) 



Note. — L is the physical domain length, 
CAi the Alfven speed. Note that the 
scaled electric field is not a dimensionless 
variable, but the inverse of a time scale. 
The obvious non-dimensionalization by tak- 
ing cEL/{\qBo) proves to be inconvenient 
because of the Ao-dependence. 



Table 2 

Time scales used in the dimensionless equations 



timcscale 


symbol 


definition 


Ohmic 


TQ 




ambipolar 


TAD 


L^/\ad 


ionization 


n 


(Ripna/Pio)~^ 


recombination 


Tree 


{RrecPio) ^ 


inductive 


Tind 


LBq/cE 


Alfven 


TA 


L/CAi 



Note. — Time scales used in the dimen- 
sionless equations (14), (15, (14). The phys- 
ical variables are defined in Table 1. 
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Table 3 

Physical parameters and time scales 



cjUciiiiiTiy 




\pj 






n„ [cni-3] 


102 


104 


10*^ 


1013 




1.3 


2.2 


2.2 


2.2 




12 


29 


29 


40 




10-4 


10-7 


10-8 


10-12 


Bo[G] 


10-5 


5 X 10-5 


2 X 10-4 


1 


T„[K] 


80 


10 


30 


500 


Csi [cm/s] 


3.0 X 10" 


6.9 X 103 


1.2 X 104 


4.1 X 104 


CAi [cm/s] 


6.3 X 10*^ 


6.4 X 107 


8.1 X 107 


1.1 X 10i° 


Csn [cm/s] 


9.2 X 104 


2.5 X 104 


4.3 X 104 


1.8 X 105 


Cyi« [cm/s] 


1.9 X lO'^ 


7.4 X 104 


3.0 X 104 


4.7 X 104 


T» [yr] 


3.4 X 10^ 


7.6 X 10° 


2.9 X 10-1 


1.6 X 10-7 




4.6 X 10-5 


2.3 X 10-8 


4.3 X 10-8 


2.9 X 10-11 


^gyro [cm] 


2.6 X 10*^ 


2.9 X 105 


1.3 X 105 


1.2 X 102 


i^D [pc] 


4.2 X 10-3 


1.6 X 10-2 


3.9 X 10-4 


2.7 X 10-7 


Lcol [pc] 


6.7 X 10-7 


2.1 X 10-9 


3.6 X 10-11 


1.2 X 10-17 



Note. — Parameters for (a) diffuse clouds, (b) dense clouds, (c) 
cores and (d) protoplanetary disks. 7 = Cp/c„ = 5/3. Molecu- 
lar weights jjL, sound speed Cg and Alfvcn velocity ca refer to the 
ions or neutrals, depending on their index i or n. The ioniza- 
tion time is given in §2.2.3, and the plasma (3 refers to the ions. 
Equations (73) and (74) define Lad and Lcoi- The ion gyroradius 

^gyro = C-JkT]I^jnH / cB . 



Table 4 
Simulation parameters 



Modol 


J 


/ } 


•AL) 




/ 


Z/m 


^34 


10-3 


104 


104 


10 


* 


4.4 X 10-11 


^44 


10-4 


104 


104 


10 


* 


7.0 X 10-13 


^32 


10-3 


102 


104 


10 


* 


4.4 X 10-13 




10-4 


102 


10* 


10 


* 


7.0 X 10-15 



Note. — Time scales (in years) and /? for the simula- 
tions. For all runs of type A we have 10^ < rn < 10i4, 
denoted by *. Sec §2.2.6 for a discussion of the physi- 
cal regimes covered. The first digit in the model name 
indicates /3, the second Tj. 
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